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We report measurements of branching fractions for charged and neutral B —> rjcK decays where 
the 77c meson is reconstructed in the K^K^'K'^,K+K-'K°,K*°K-Tr+ and pp decay channels. The 
neutral B° channel is a CP eigenstate and can be used to measure the CP violation parameter 
sin 2^1. We also report the first observation of the B" rjcK*" mode. The results are based on an 
analysis of 29.1 fb~^ of data collected by the Belle detector at KEKB. 



PACS numbers: 13.20.He, 13.25.Hw, 13.60.Rj 
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The decay mode B — > rjcK proceeds by a spectator 

b — > CCS transition with internal PF-cmission as in the 
CP eigenstate B° J/tpKg. The neutral decay mode 
— > rjcKg has therefore been used to measure the CP 
violation parameter sin 2(^1 [1], [2], [3]. Measurements of 
branching fractions for B rjcK^*^ decay modes are also 
useful in the study of the dynamics of hadronic B decay 
[4]. However, in contrast to i?" J/tpKg, the rjc meson 
must be reconstructed from hadronic decays rather than 
from a leptonic final state with relatively low combinato- 
rial background. In this paper, we report new measure- 
ments oi B ^ rjcK branching fractions with the 77c me- 
son reconstructed in the K^K+tt' , K' K+n^ , K*^K-'!r+ 
and pp channels [5] . These signals are large enough to be 
used to determine the mass and width of the rjc meson. 
We also report the first observation of the related decay 
mode B° ^ r?ci^*°- When K*'^ K^n'^, this decay 
mode is also a CP eigenstate. 

We use a 29.1 fb~^ data sample, which contains 31.3 
million produced BB pairs, collected with the Belle de- 
tector at the KEKB asymmetric-energy e+e^ (3.5 on 
8 GeV) collider [6]. KEKB operates at the T{4S) res- 
onance {y/s = 10.58 GeV) with a peak luminosity that 
now exceeds 7 x 10^'^ cm^^s^^. The Belle detector is 
a large-solid-angle magnetic spectrometer that consists 
of a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), a mosaic of aerogel thresh- 
old Cerenkov counters (ACC), time-of- flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located 
outside of the coil is instrumented to identify and 
muons (KLM). The detector is described in detail else- 
where [7]. 

We select well measiired charged tracks with impact 
parameters with respect to the interaction point of less 
than 0.5 cm in the radial direction and less than 3 cm 
in the beam direction (z). These tracks are required to 
have pt > 50 MeV/c. 

Particle identification likelihoods for the pion and kaon 
particle hypotheses are calculated by combining infor- 
mation from the TOF and ACC systems with dE/dx 
measurements in the CDC. To identify kaons (pions), 
we apply a mode-dependent requirement on the kaon 
(pion) Hkelihood ratio, Lk / [Ltt + Lk) {Lt^ / {Lt^ + Lk))- 
The requirement LxliL-n -I- Lk) > 0.5 is used for the 
rjc KgK^TT^ mode. For other modes, we require 
Lk/{Lt, + Lk) > 0.6, which is 88% efficient for kaons 
with a 8.5% misidentification rate for pions. For the 
r]c K~^K~Tr° mode, the kaon likelihood ratio is re- 
quired to be greater than 0.8 for those charged kaons that 
come directly from the B rather than from the rjc can- 
didate. In addition, we remove all kaon candidates that 
are consistent with being either protons or electrons. 



Protons and antiprotons are identified using all particle 

ID systems and are required to have proton likelihood 
ratios {Lp/{Lp + Lk) and Lp/{Lp + L^)) greater than 
0.4. Proton candidates that are electron-like according 
to the information recorded by the CsI(Tl) calorimeter 
are vetoed. This selection is 95% efficient for protons 
with a 12% kaon misidentification rate. 

We select Ks candidates from tt+tt" candidates that 
lie within the mass window 0.482 GeV/c^ < M(7r+7r~) < 
0.514 GeV/c^ (±4a-). The fiight length of the Ks is re- 
quired to be greater than 0.2 cm. The difference in the 
angle, in the x — y plane, between a vector from the beam 
spot to the Ks vertex and the Ks fiight direction is re- 
quired to satisfy A(f> < 0.1 rad. 

K*^ candidates are reconstructed in the K^n^ mode. 
For r]c — > K*^K~Tr~^ , we require the X+tt" invariant 
mass to be between 0.817 and 0.967 GeV/c^. For the 
B^ —* rjcK*'-' mode, the K*'^ mass must lie in the range 
between 0.801 and 0.991 GeV/c^. 

Neutral pion candidates are selected from pairs of ECL 
clusters with invariant mass within ±16 MeV of the nom- 
inal 7r° mass and momenta above 350 MeV/c. The pho- 
tons must have energy above 50 MeV if they lie in the 
barrel region of the calorimeter and above 200 MeV if 
they are detected in the endcap. 

To reconstruct signal candidates in the _B+ — > rjcK'^ 
and B^ — s- rjcK'^ modes, we form combinations of charged 
or neutral kaons and t]c candidates. The r]c is recon- 
structed in the K^K+tt' , K" K+tt° , K*° K' tt+ and pp 
decay modes. The ric candidate is required to have in- 
variant mass in the range 2.920 < Af^^ < 3.035 GeV/c^ 
for the K-K+TT^ mode and 2.935 < < 3.035 GeV/c^ 
for all other modes. The charged daughters of the rjc are 
required to satisfy a vertex constrained fit with a mode- 
dependent requirement [8]. 

To isolate the signal, we form the beam-energy con- 
strained mass, Mbc = \/ E^eam - Brecon: and energy dif- 
ference AE = i^recon — -Ebeam in the T(4S') Center of mass 
frame. Here -Bbcam, £'rocon and Prccon are the beam en- 
ergy, the reconstructed energy and the reconstructed mo- 
mentum of the signal candidate, respectively. The signal 
region for AE in all modes except for rjc — > KgK~TT~^ 
is ±2.5(7, where a is the mode-dependent resolution 
and ranges from ±25 MeV for r]c pp to the range 
-55, +45 MeV for r/^ ^ K-K+tt^. In the low back- 
ground rjc — ^ K^K^t:^ mode, the region is extended 
to ±35 MeV, (±3.5a-). The signal region for Mbc is 
5.270 GeV/c2 < Mbc < 5.290 GeV/c^. The resolution 
in beam-energy constrained mass is 2.8 MeV/c^ and is 
dominated by the beam energy spread of KEKB. 

Several event topology variables provide discrimina- 
tion between the large continuum (e+e^ qq, where 
q = u,d, s, c) background, which tends to be coUimated 
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along the original quark direction, and more spherical 
BB events. We first remove events with R2 > 0.5, where 
i?2 is the normalized second Fox- Wolfram moment. We 
form a likelihood ratio using two variables. Six modi- 
fied Fox- Wolfram moments and the cosine of the thrust 
angle are combined into a Fisher discriminant [9]. For 
signal MC and continuum data, wc then form probabil- 
ity density functions for this Fisher discriminant, and 
the cosine of the B decay angle with respect to the z 
axis (cos6'b). The signal (background) probability den- 
sity functions are multiplied together to form a signal 
(background) likelihood jCs {^bg)- A mode-dependent 
likehhood ratio requirement Cs/{Cs + C-bg) is then im- 
posed. 

Using a sample of 57 million BB Monte Carlo events 

with a model of 6 ^ c decays, we investigate backgrounds 
from other B decay modes. In the rjc pp mode, no 
such backgrounds are found. In other modes, some back- 
ground is observed but it can be removed by application 
of mode-dependent vetos on invariant mass combinations 
that are consistent with the D, Dg, Xci, J/'^, i'i'^S) 
or r]c{2S) masses. For example, in the B~^ rjcK^, 
rjc K^K^TT^ mode, we find there is background 
from the decay chain B+ DOp+, D° K-K+, 
tt+tt". This background is removed by requiring 
that the invariant mass be inconsistent with the 

mass. 

We fit the Mbc distribution to the sum of a signal Gaus- 
sian and a background function that behaves like phase 
space near the kinematic boundary [10]. The width of the 
Gaussian is fixed from MC simulation while the mean 
is determined from i?+ Z)°7r+ data. The shape pa- 
rameter of the background function is determined from 
Ai? sideband data. The signal yield was determined by 
fits to the individual Mbc distributions for each mode. 
The yields and significances [11] for these fits are given 
in Table I. Significant signals arc observed in all decay 
modes except for B^ VcK^ , Vc ^ K*^K^'k^ . For this 
mode, we calculate an upper limit based on the num- 
ber of events observed in the Mbc signal region (4) and 
the expected number of background events (2) based on 
the fit. We use the Feldman-Cousins procedure [12], and 
reduce the efiiciency by one sigma of the systematic er- 
ror in the calculation [13]. The detection efficiencies for 
all modes were determined from a GEANT based Monte 
Carlo simulation. 

For illustration, in Fig. 1, we show the beam-energy 
constrained mass and Ai? distributions for the signal can- 
didates in all the decay modes except for B^ VcKg, 
r]c — » K*°K~'K~^ . In the Mbc distribution, we observe a 
signal of 195 ± 17 events. 

As a consistency check, we also determine the yield 
from a fit to the AE distribution with a double Gaussian 
for signal and a linear background function with slope 
determined from the Mbc sideband. The results of these 
fits are also given in Table I. The fit to the AE distri- 




M;„(GeV/c^) AE(GeV) 
FIG. 1. (a) Mbc and (b) AE distributions for B -» r]cK 
candidates in all the decay modes. 

bution for all modes combined gives an integrated yield 

of 188 ± 17 events in the signal region. In the fits to the 
AE distribution, the region with AE < —120 MeV is ex- 
cluded to avoid contributions from modes with additional 
particles such as B ^ VcK* . 

After removing the requirements on rjc invariant mass, 
we also verify that the signal yield for r]c candidates in 
the Mbc, AE signal region is consistent with the re- 
sult used for the branching fraction determination. The 
ric invariant mass distribution for signal candidates is 
shown in Fig. 2. Fitting to a Breit-Wigner convolved 
with the resolution determined from MC, we find an 
intrinsic width r{r]c) = 29 ± 8 ± 6 MeV and a mass 
M{r]c) = 2979.6±2.3±1.6MeV. The systematic errors in 
the width and mass measurements include the effects of 
varying the background shape, the small difference be- 
tween data and MC detector resolutions, and possible 
binning effects. The results are consistent with world av- 
erages and comparable in precision to the best individual 
measurements [14]. The yield is 182 ±25 events. We also 
observe a clear signal of 66 ± 18 events from B J/tpK, 
where the J/tp is reconstructed in hadronic decay modes, 
that has a J/iJj mass and yield that are consistent with 
values obtained for the J/t/j (.'^(r decay mode. 
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M(ii,)(GeV/c^) 

FIG. 2. Candidate Mijjc) invariant mass distribution for 

events in the Mbc and AE signal region. Signals at the 77c 
and J/tp masses from B — > ricK and B —> J/ipK decays are 
visible. 

The contributions to the systematic error include the 
uncertainties due to the tracking efficiency (2% per 
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track), particle identification cfRcicncy (4-14%, depend- 
ing on the mode) and the modehng of the hkelihood 
ratio requirement (2%). The error in kaon identifica- 
tion efficiency is obtained from kincmatically selected 
£)*+ £)°7r+, K~Tr^ in the data while the error in 

proton/ antiproton identification is determined using A/ A 
samples. The systematic error due to the modeling of 
the likelihood ratio cut is determined using i?+ l)°7r+ 
events reconstructed in data. The systematic error in 
the yields of the Mbc fit were determined by varying the 
mean and a of the signal and the shape parameters of 
the background. To account for the possibility of back- 
ground from non-resonant modes that may contribute 
to the Mbc distribution, we include the yields observed 
in the r]c mass sideband (8-14% of the signal depending 
on the mode) as an asymmetric systematic error. The 
sources of systematic error are combined in quadrature 
to obtain the total systematic error, which is given in 
Table II. 



and in some cases there arc conflicting measurements, 
we quote B branching fractions for the rjc KgK~n~^ 
and r]c — > K~K~^'k^ modes only. The r]c — > KgK~TT~^ 
mode is the most precisely and reliably measured mode 
and the branching fraction for the rjc K~ K~^t^^ mode 
is related by isospin. We use B{ric — > K^K~-k^) = 
1/3 X (0.055 ± 0.017) where 1/3 is the product of the 
appropriate Clebsch-Gordon coefficient and intermedi- 
ate branching fraction. We assume that the exper- 
imental systematic errors in the 77c KlK-TT+ and 
r]c K~K~^TT° modes are uncorrelated. We assume 
equal production of B^B~ and B^B^ pairs and do not 
include an additional systematic error for the uncertainty 
in this assumption. We find 



B{B+ -> ricK+) = (1.25 ± 0.141^;}^ ± 0.38) x 10 



and 



B(B" ^ r?c^") = (1-23 ± 0.23tg;i^ ± 0.38) x 10 



TABLE 1. Signal yields from Mbc and fits, statistical 
significances and MC reconstruction efficiencies. Errors are 

statistical only. 



B+ ^ 'q^K+ 


AE Yield 


Mbc Yield 


Signif.(Mbc) 


e{%) 


rjc 




74.8 ± 10.4 


81.6 ±10.3 


12.2cr 


16.4 


r]c 




26.5 ± 7.8 


31.8 ±7.0 


6.3a 


8.8 


Vc 


PP 


16.3 ±4.2 


17.7 ± 4.4 


7.5cT 


34.0 




K*°K-n+ 


22.0 ± 5.8 


20.8 ± 5.4 


5.8cT 


5.1 



B°- 




AE Yield 


Mbc Yield 


Signif.(Mbc) 


e{%) 


Vc 




KlK-iT+ 


19.6 ± 5.4 


23.0 ± 5.4 


6.8cr 


15.5 


Vc 






19.9 ±5.7 


17.1 ±5.1 


4.7a 


9.5 


Vc 




PP 


7.0 ±3.0 


6.8 ±2.6 


5.0a 


34.9 


Vc 






0.2 ±1.7 


2.2 ±1.8 


1.6a 


3.75 



TABLE 11. Product branching fractions for B 

cay modes (10~®). 



rjcK de- 



B{B+ 
B{B+ 
B{B+ 
B{B+ 



VcK+) X S(77c 
??cA'+) X S(r?c 
'ncK+) X B{ric 
r,,K+) X S(»?c 



■7r+) 



B(B° -> ri^K") X B(r;c 
B(S" ^ r/cAT") X S(r;c 
B{B" ^ r/cit") X S(r?c 
B{B° ricK°) X B{r]c 



K+K 

-^PP) 

■ K°sK-n+ 

■ PP) 



^-3.8/ 
1.1^ 



(23.2 ± 2.t_ 
(11.4±2.5liij 

(19.3 ±5.01^3:^) 

{20.1±4.Ttl;l) 
(16.6±5.0ii-^) 
(1.79±0.68l«:^») 
(8.1 ±6.6 ±1.4) 
< 29 at 90% C.L. 



The product branching fractions are given in Table 
II for all modes in which signals are observed. Since 
many of the r]c branching fractions are poorly determined 



The first error is statistical, the second error is system- 
atic and the third error is due to the uncertainty in the 
r]c branching fraction scale. When the rjc branching frac- 
tions for the other modes are better determined, absolute 
B branching fractions for these modes can be extracted 
from our results. 

In the B° VcK*°, K*^ -» AT+tt" channel, the r?c is 
reconstructed in the K^K^tt^ mode. Since this mode is 
a pseudoscalar to pseudoscalar-vector decay, by angular 
momentum conservation, the cosine of the K* helicity 
angle {cos Oh) follows a cos^Oh distribution; we select 
events with | cos9h\ > 0.4. We also investigate BB back- 
ground and find that this background can be removed 
by applying vetos to events with combinations that are 
consistent with J/i/i K^Kir, J/ip K+R- 
V'(25) ^ K^sKtt, 7^23) ^ KlKn, xa K^K^, 
rj, K+R-TT+TT-, Df KlK+ and D, -f K-R+n. 
The detection efficiency for these selection requirements 
is (7.95 ±0.12)%. 

After applying these requirements, a fit to the Mbc 
distribution yields a signal of 33.7 ± 6.7 events for B^ —^ 
rjcK*'^ with a statistical significance of 7.7(t. [11] The Mbc 
distribution is shown in Fig. 3(a). The yields from the 
AE fit, shown in Fig. 3(b) (30±7 events), the rjc invariant 
mass distribution (24 ±7 events) and the /T+tt" invariant 
mass distribution (27 ± 8 events) are consistent with the 
yield from the A'/bc fit. 

To evaluate the contribution from non-resonant B'^ —> 
KgK'^TT~K* as well as the remaining BB backgrounds 
that peak in the Mbc distribution, we select events in 
the rjc sideband [15] and repeat the Mbc fit- We find 
no significant signal. By using the ratio of the yields in 
the r/c signal and sideband regions determined from MC, 
we estimate the contributions of such backgrounds to be 
3.9±4.2 events, consistent with zero. We use the K* side- 
band [15] to estimate the non-resonant B^ rjcKn decay 
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component and obtain —0.6 ± 3.3 events. These possi- 
ble background contributions are not subtracted in the 
branching fraction measurement, but instead are treated 
as systematic uncertainties. We find 

B{B° VcK*°) = (1.62 ± 0.32t° ± 0.50) x 10"^. 

To take into account the possibihty of BB background, 

wc conservatively include asymmetric systematic errors 
from the results of the fits to the ryc(— 12%) and K* 
(—8%) sidebands. Other sources of systematic error are 
the uncertainties in the track reconstruction efficiency 
(±2% per track), the parameters in the Mbc fit (± 7%), 
particle identification (±6%) and the number of B^B^ 
events. 

From the results for the branching fractions for 5° 
rjcK^ and B^ — > r]cK*^ determined above, we can deter- 
mine the ratio i?,,,^ = B{B" r]eK*")/B{B" VcK°). 
The uncertainty from the rjc branching fraction scale can- 
cels in the ratio. We find 

Rr,, = 1.33 ±0.36^^;^;^. 

Our result can be compared to calculations of this ratio in 
models based on factorization and is consistent with the 
range 1.02-2.57 predicted by Gourdin, Keum and Pham. 
[16]. 
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FIG. 3. (a) Mbc and (b) AE distributions for B rjcK*° 
candidates. 

The branching fractions reported in this paper for 
5+ VcK^ and B° — » rjcK^ using rjc KKw decays 
are more precise than previous results [17]. The result for 
the _B+ rjcK^ branching fraction is somewhat higher 
than the CLEO measurement, while the B^ rjcK^ re- 
sult is consistent. Several additional r]c modes including 
T]c ^ PP, Vc ^ K~K'^iT^ and rjc K*'^K~7r^ have 
been used and increase the fraction of rjc decays that can 
be reconstructed for CP violation measurements. With 
the large samples of _B ^ rjcK decays now available, 
we are able to determine the mass and width of the r]c 
meson. We find M{r]c) = 2979.6 ± 2.3 ± 1.6 MeV and 
r{ric) = 29 ± 8 ± 6 MeV. In addition, we report the first 
observation of the 5° rjcK*^ decay, which is an CP 
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